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Abstract 

J> ■ We present a general framework to study the effect of killing sources 

. on moving particles, trafficking inside biological cells. We are merely 

| concerned with the case of spine-dendrite communication, where the 

CN ' number of calcium ions, modeled as random particles is regulated across 

the spine microstructure by pumps, which play the killing role. In par- 
[ ticular, we study here the survival probability of ions in such environ- 

' ment and we present a general theory to compute the ratio of the number 

Qh 1 of absorbed particles at specific location to the number of killed particles 

^ fii during their sojourn inside a domain. In the case of a dendritic spine, 

the ratio is computed in terms of the survival probability of a stochastic 
trajectory in a one dimensional approximation. We show that the ratio 
depends on the distribution of killing sources. The biological conclusion 
. ^ | follows: changing the position of the pumps is enough to regulate the 

^ ■ calcium ions and thus the spine-dendrite communication. 
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1 Introduction 



The post-synaptic part of a synapse is usually a dendritic spine, a microstruc- 
ture located on a dendrite of a neuron (see figure [TJ) |IT|, [T2], [TU] . The spine 
geometry consists of a nearly spherical head connected to the dendrite by a 
narrow cylindrical neck. Calcium ions enter the spine head through gluta- 
mate gated channels following the release of glutamate neurotransmitters by 
the pre-synaptic terminal. The communication between a dendritic spine and 
the dendrite depends on the ability of the calcium ions to pass through the 
cylindrical neck to the dendrite. When ions enter the neck, they diffuse and 
either reach the dendrite or are extruded on their way to the dendrite by pump 
proteins located on the lateral surface of the neck f3J ■ 

The number of calcium ions that arrive at the dendrite and the calcium 
contents of the spine are regulated by the geometry of the neck and by the 
contents of the spine. The contents include organelles, such as the endoplasmic 
reticulum, calcium buffer proteins such as calmodulin, calcium stores, actin- 
myosin proteins, and pumps on the spine membrane. In this paper, we focus 
on the role of the spine neck in spine-dendrite communication, which is an area 
of intense experimental research (see, e.g., [TI], |12j). 

We adopt a simplified one-dimensional model of the diffusive motion of 
calcium ions in the neck, in which the termination of ionic trajectories by 
pumps is described as "killing", while termination in the dendrite is described 
as "absorption". A killing measure is the probability per unit time and unit 
length to terminate a trajectory at a given point at a given time. Thus an 
ion can pass through a killing site many times without being terminated. In 
contrast, an absorbing boundary terminates the trajectory with probability 
1 the first time the trajectory gets there. Thus we distinguish between two 
random times on a trajectory, the time to be killed, denoted T, and the time 
to be absorbed, denoted r. 

We need to find the probability Pr {r > T | y} of an ion getting killed 
(pumped out) in the neck before it is absorbed at the boundary (the den- 
drite), given that it started at a point y in the neck. The ratio 

Pr{r <T\y} 
00 Pi{r>T\y} 

is the fraction of absorbed to killed (pumped) particles. We also need to 
calculate E [T | r > T, y] , the mean time to be killed, given that the particle is 
killed, as well as E [r | T > r, y] , the mean time to absorption, given that the 
particle is absorbed. 

An application of our model in neurobiology concerns calcium regulation 
in the dendritic spine and in the dendrite. In dendrites of neurons, ions are 
constantly exchanged between compartments and when the concentration of 
calcium ions in the dendritic shaft rises above a threshold value, some specific 
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cascades of chemical reactions are initiated, that can lead to a new physio- 
logical stage, where the synaptic properties are modified. For example, the 
biophysical properties of the synapses or the number of channel receptors can 
be irreversibly changed PH1IE1IH]- The process that consists of changing the 
synaptic properties is known as synaptic plasticity. Today, the mechanisms of 
induction of synaptic changes are still unclear, but it has been demonstrated 
recently ^H] that the induction process can be affected by the dynamics of the 
spine-dendrite coupling. The communication between a dendritic spine and 
the dendrite depends on the ability of the ions to pass through the cylindrical 
neck of the spine (see figure Q). The measure of this ability is the parameter 
Roo. When ions leave the spine head and enter the neck, they diffuse and 
either reach the dendrite (with probability Pr {T > r | y}), or, as mentioned 
above, are extruded by pump proteins on their way to the dendrite [T3J[TT]. 

In a simplified homogenized model proposed in [16] . the number of ions 
filtered by the neck has been estimated and compared with experimental re- 
sults. This number depends on the distribution of pumps along the neck and 
on the efficiency of the pumping process. The precise comparison with exper- 
imental data in ^B] made it possible to predict that changing the length of 
the spine neck (which occurs under specific conditions, see for example |I2j) 
is sufficient to regulate precisely the number of ions arriving at the dendrite. 
Spine-dendrite calcium signaling (OIH1) and its regulation through specific 
microstructures, such as the spine neck, is crucial for the induction of synaptic 
plasticity, which underlies learning and memory. 

The mean time E[t \ y] an ion spends inside the neck can be written as 

E[t\y] = E[t\T<T,y]Pi{r<T\y} + E[t\T<T,y]Pi{T<r\y} 

= E[t\t <T, y] Pr{r <T\y} + E[T \T <r,y] Pr{T < r | y}. 

The rate 1/E[t\y] is the total probability flux out of the neck. This is a 
measurable quantity that can be used to prove that ions diffusing into the 
dendrite originate in the spine head. Indeed, calcium that enters the spine 
head through the glutamate gated channels at the top of the spine head takes 
much longer to reach the dendrite than calcium that enters through voltage 
gated channels. This is due to the much faster propagation of the membrane 
depolarization than movement by diffusion. 

In a biological context the final distribution of particles between absorption 
and killing indicates the future changes in the steady properties of the synapse. 
This is a general principle in cell biology regulation. It is fundamental for 
the homeostasis of a living cell to regulate the number of proteins or small 
molecules it contains and to maintain this number constant in the absence of 
external input. This is for example achieved through an equilibrium between 
synthesis and hydrolysis mechanisms. At a molecular level, when molecules 
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Dendrite 



Dendrite 




Figure 1: Extrusion of an ion from the spine neck. A dendritic spine is 
a microstructure located on the dendrite of neurons, consisting of a round head 
connected to the dendrite through a cylindrical neck. Its function is still unclear. 
After ions enter through the head, either they are pumped out (right figure) or 
they reach the dendrite (left picture). The number of ions reaching the dendrite is 
regulated by the number and the distribution of pumps and the length of the spine 
neck. The neck length changes dynamically and this is induced by previous calcium 
ions. 

reach the active sites of free enzymes by a Brownian random walk, either 
the molecules are hydrolyzed or nothing happens (see fH] for a stochastic 
description) and after some time, the molecules are absorbed or enter different 
organelles. This is what happens in signal transduction, as in synapses of 
neurons or in sensor cells. In some cases, the stability and the function of the 
cell depends on the efficiency of such dynamical processes. In addition, the 
geometry of the cell participates in the regulation of the number of particles, 
such as ions, that reach specific locations. 

In the present work, we are interested in estimating the probability that 
an ion survives in a medium containing many pumps. We compute the prob- 
ability to arrive at a specific location before being killed (see figure EJ) as a 
function of structure and pump distribution. In the case of a dendritic spine, 
we assume that the cylindrical neck can be approximated by a one dimen- 
sional interval, and the computations are given in a one-dimensional model. 
The one-dimensional approximation is valid when the radius of the spine neck 
is sufficiently small, otherwise, the small pumps cannot affect the normal dif- 
fusion process (see |2|). We will see that various pump distributions affect the 
concentration of ions in the neck; we compare a uniform distribution along the 
spine neck, modeled as a constant killing rate, to an accumulation of pumps in 
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"hot spots" at some specific locations, (for example at the base of the dendritic 
spine). In either case, we estimate the flux of ions into the dendrite. 

The reduced one-dimensional model of Brownian motion with killing and 
absorption is investigated in various types of killing sets. It is of interest to de- 
termine the influence of spatial distribution of the killing measure on the global 
survival probability of the population. Absorption and killing are expressed 
differently in the Fokker-Planck equation (FPE) for the transition probability 
density function (pdf) of the Brownian motion. While total absorption at the 
boundary is expressed as a homogeneous Dirichlet boundary condition, killing 
appears as a reaction term in the FPE pQ. 

Our main results are general expressions for the probabilities, ratio, and 
mean times in general, and in particular, we give explicit expressions as func- 
tions of the geometry and distribution of killing sites in the one-dimensional 
model. We also provide a biological interpretation of the results. 

2 Killing measure and the survival probability 

We consider a Brownian motion (particle) in a cylinder, whose lateral bound- 
ary contains many small absorbing hole, one base is reflecting and the other 
absorbing. This model can be approximated [2] by a one-dimensional Brown- 
ian motion on an interval with one reflecting and one absorbing endpoints, and 
a killing measure inside the interval. The strength of the killing measure is re- 
lated to the absorption flux of the three-dimensional Brownian motion through 
the small holes on the boundary of the cylinder. The killing measure k(x, t) is 
the probability per unit time and unit length that the Brownian trajectory is 
terminated at point x and time t |XJ. 

The survival probability and the pdf of the surviving trajectories can be 
derived from the Wiener measure [3]. Indeed, for a Brownian trajectory X(t) 
and the random time at which it is terminated, r, we denote the (defective) 
probability density of finding a trajectory at point x, given that it starts at y, 
by 

u(x, 1 1 y) dx = Pr {x(t) Gil dx, r > t \ x(0) = y} . 

To derive the joint density of x{r) and r, we can formulate the problem in 
terms of the Wiener integral with a killing measure. The Wiener density per 
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unit time of being killed in the time interval [t, t + At] at a point x^ = x is 
Pt{x(t) =x, t = t\x(0) =y} = (2.1) 



lim — k (x N , t N ) At 

N^oo /\ t 



N 



N 

IJ ( exp 

3=1 



[Xj Xj-lj 



2At 

k (x, t) u(x,t\y) , 



[1 — k (xj-i, tj-i) At] dxj-i 



where 



t 



At = -, t 3 = 3 At, 

and u(x,t\ y) is the solution of the initial value problem 

Ut = u xx — k (x, t) u, for x G R, t > 

u (x, 0) = 5 (x — y) . 
In the case that k (x,t) = V and the diffusion coefficient is D, we have 

du(x, t\y) 3 2 u(xA\y) Tr . « „ m 

= D \ '„ 1 yj - V u(x, t\y), for x e R, t > 



(2.2) 



dt 



dx 2 



(2.3) 



u(x,Q\y) = 8(y-x). 
The solution is given by 

1 

u(x,t\ y) 



2VM 



exp <^ -V t 



(x - yf 

ADt 



(2.4) 



The effect of absorption is expressed through many different features of the 
Wiener integral. First, the probability per unit time of being killed (absorbed) 
inside the interval [a, b] at time t is 



Pr {x(t) E [a,b], t = t \ x(0) = y} 



k (x, t) u(x,t\ y) dx, 



while the probability of being killed in the interval before time t is 



Pr {x(t) E [a,b], t < t \ x(0) = y} 



t f b 



Ja 



k(x, t)u(x, t\y) dx dt. 



The probability of ever being killed in the interval is 

-oo pb 



poo pb 

Pr {x(t) E [a, b] | x(0) = y} = k{x,t)u{x,t\y) dxdt, 

JO J a 

and the density of ever being killed at x is therefore 

POO 

Pr{x(r) = x\ x(0) = y} = / k{x, t)u (x, t \ y) dt. (2.5) 



o 



The survival probability is the probability that the trajectory still exists at 
time t, that is, 

S(t) = Pr{r > t\ x(0) = y} = / u(x,t\y)dx. 

m 

For the case k(x,t) = V eq. ()2.4|) gives 

Pr{r >t\x(0) =y}= / u(x, 1 1 y) dx = e~ Vot . (2.6) 

JR 

This is exactly the rate at which particles disappear from the medium. The 
rate is exponential, so that out of No initial independent Brownian particles 
in M the expected number of particles that have disappeared by time t is 
N (l — e~ Vot ). The probability of being killed at point x, given by eq. ()2.5|) . is 



1 v„ v a . 

= -\ — exp < — \ — \x — y\ 

2 V D F | VD 1 y| 

We assume henceforward that the killing measure is time independent. 



(2.7) 



3 Absorption versus killing 

We consider now a particle diffusing in a domain Q C lR n with a killing measure 
k(x) and an absorbing part dQ a C dQ of the boundary dQ. Thus the trajectory 
of the particle can terminate in two ways, it can either be killed inside Q 
or absorbed in dQ a . The difference between the killing and the absorbing 
processes is that while the trajectory has a finite probability of not being 
terminated at points x where k(x) > 0, it is terminated with probability 1 
the first time it hits dQ a . Thus the trajectory may traverse many times killing 
regions, where k(x) > 0, but it cannot emerge from the absorbing part of the 
boundary. 
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3.1 Definition and basic equations 

We define two random termination times denned on the trajectories of the 
diffusion process: the time to killing, denoted T, and the time to absorption 
in dfl a , denoted r, which is the first passage time to dfl a . We calculate below 
the probability Pr{T < r | y}, and the conditional distribution Pr{r < t\r < 
T,y}. 

We consider the trajectories of the stochastic differential equation 

dx = a(x)dt + B(x)dw(t) for x(t) G 0, (3.1) 

where a(x) is a smooth drift vector, B(x) is a smooth diffusion matrix, and 
w(t) is a vector of independent standard Brownian motions pQ. We assume 
that a killing measure k(x) > is defined in O and k(x) > on a set of 
positive measure. 

The transition probability function of x(t) satisfies the Fokker-Planck equa- 
tion 

Op(x t yj 

= £p(x, 1 1 y) - k(x)p(x, 1 1 y) for x, y G O, (3.2) 

where the forward operator C is defined by 

r„(„ + \..\ _ d 2 a^(x)p(x, t \ y) ^ da l {x)p{x,t\y) 

L P {X,t\y)-2^ g-g- 2^ Q? > W 

t,j=l 1=1 



and 



a(x) = l -B{x)B T {x). 



The forward operator L can also be written in the divergence form 

Cp(x,t\y) = -V- J(x,t\y), (3.4) 
where the components of the flux density vector J(x, 1 1 y) are defined as 

J'(x,t\y) = -J2 d "' J(x) ^ My) + a>(x) P (x, 1 1 ,). (3.5) 

3=1 X 

The initial and boundary conditions for the Fokker-Planck equation ()3.2|) are 
p(x,0\y) = 5(x — y) for x,yeil (3.6) 

p(x,t\y) = for t > 0, x G 90, y G tt a (3.7) 

J(x,t\y)-v(x) = for t > 0, x G 90 - <90 a , y G O. (3.8) 
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The transition pdf p(x, t\y) is actually the joint pdf 

p(x, t\y)dx = Px{x(t) E x + dx, T > t, r > t\y}, (3.9) 

that is, p(x, 1 1 y) is the probability density that the trajectory survived to time 
t, i.e., was neither killed nor absorbed in dfl a , and is located at x. 
We begin by showing that 

poo P 

Pr{T < r\y} = / / k(x)p(x,t \ y) dxdt (3.10) 
Jo Jn 

by two different derivations. First, assume that the entire boundary is ab- 
sorbing, that is, dQ a = dQ. Then the probability density of surviving up to 
time t and being killed at time t at point x can be represented by the limit as 
iV -> oo of 

Pr $x N (ti )N ) e O, x N (t 2jN ) efl,..., x N (t) = x,t <T <t + At | x(0) = t/j 




>n Jn Jn 



N 

n 



\ ^(27rAt) n det ( r(x)(t j _ hN )) 



1 rp 

x exp 1 ~ 2At ~ x (Ji-^) ~ a(a;(tj_i ) jv))At] cr~ (x(tj^ ltN )) 



x 



fe(aj)At, (3.11) 



[y j - x(tj- ljN ) - o(aj(tj-_i,jv))At] | [1 - fc(aj(fy,jv)At] 
where 

At = ^, *i,iv = J At, 

and 

as(*o,Jv) = 2/ 

in the product. The limit is the Wiener integral defined by the stochastic 
differential equation (|3.1|) . with the killing measure k(x) and the absorbing 
boundary condition [T2|. In the limit iV — > oo the integral (J3.ll)) converges 
to the solution of the Fokker-Planck equation ()3.2j) in Q with the initial and 
boundary conditions ()3.6|) and ([3.7)1 . Integrating over f2 with respect to a? and 
from to oo with respect to t, we obtain, in view of tj3.9j) . the representation 

A second derivation begins with the integration of the Fokker-Planck equa- 
tion (13~21) . 

/'OO /» /*OG /» 

1= / * J(x,t\y)'u(x)dS x dt+ / / k(x)p(x,t\y) dxdt. (3.12) 



9 



We write 



J{t | y) = I J{x, t\y)- u{x) dS x (3.13) 
Jan 

and note that this is the absorption probability current on dQ. Therefore, in 

/oo 
J(t I y) dt is the total proba- 
o 

bility that has ever been absorbed at the boundary dVt a . This is the probability 
of trajectories that have not been killed before reaching dQ a . Writing eq. (|3.12|) 



as 



/*oo roc r 

/ J(t | y) dt = 1 — / / k(x)p(x, 1 1 y) dx 
Jo Jo Jq 



dt, 



we obtain f)3.10j) . 

The probability distribution function of T for trajectories that have not 
been absorbed in dQ a is found by integrating the Fokker-Planck equation with 
respect to x over Q and with respect to t from to t. It is given by 



Pr{T < 1 1 r > T, y} 



Pr{T < t,r > T\ y} 
Pr{r > T\y} 



o Jn 



k(x)p(x, s\y) dx ds 



JQ, 



k(x)p(x, s\y) dx ds 



(3.14) 



Hence 



OO POO 



k(x)p(x, s\y) dx ds dt 

E[V\ I < t. y] = J - % J r n . (3.15) 

k(x)p(x, s\y)dx ds 



o Jo, 



Equivalent ly, 



/oo p 
s I k(x)p(x, s \y) dx ds dt 

r. [ > | / < 7 . „ ° r J n f , (3.16) 

k(x)p(x, s | y) dx ds 



o Jq 

which can be expressed in terms of the Laplace transform 



/oo 
p(x,s | y)e~ qs ds 
o 
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as 



J 


k(x)p'(x, q 

n 


y) clx 


J 


1 k(x)p(x,q 
u 


y)dx 



E[T \T <r,y) = 

yk(x)p(x, q | y)dx 
u 

(3.17) 



d_ 
dq 



In | J k(x)p(x, q \ y)dx 



q=0 



The conditional distribution of the first passage time to the boundary of tra- 
jectories, given they are absorbed, is 



J 


fj{s 




y) ds 


'-.) 


poo 
J 


k(x)p(x, s 

n 


y) dx ds 



Pr{r < 1 1 T > r, y} = . (3.18) 



Thus the mean time to absorption in dQ a of trajectories that are absorbed is 
given by [TH] 

POO 

E[t \ T > t, y] = / Pr{r > t \T > r,y}dt 
Jo 



(3.19) 



J 


i>oo 

sJ(s 




y) ds 




poo 

J 


/ k(x)p(x,s 
n 


y) dxds 



The survival probability is given by 

S{t\y) = / p(x,t\y)dx. (3.20) 

3.2 An application: a "hot spot" killing in a finite in- 
terval 

We provide here first an explicit estimation of the survival probability when 
the killing measure is a Dirac killing at a single point in a finite interval and 
second, we estimate the conditional mean first passage time to exit before 
being killed. 



3.2.1 Explicit decay of the survival probability 

To compare the survival probability of Brownian motion with and without a 
Dirac killing at a point X\ in the interval [0, 7r] with absorbing boundaries, we 
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consider the solution of the boundary value problem 

du(x,t\xi) B 2 u(xA\xi) rrr/ . , . ^ 

— v 1- 1 u = D v ^' ' -V5{x - x l )u{x ) t\x l ) onR (3.21) 

■u(x, I Xi) = <5(x — y). 

w(0, = m(7T, t I Xi) = 0, 

and we denote by G the Green function of the free particle problem, where 
V = 0, then 

2 °° 

G(x, £ | y) = — y sin nx sin m/e~ n t . 

n=l 

Therefore the survival probability of Brownian motion in the interval is 

Sb(t I ») = f G( X , t\y)d X =*± sin(2n-l )!/ e _ (2 „_ 1)a , 
Jo t Jrf 2n. - 1 

Using the Laplace transform, the solution u(x,t\ y) of ()3.21|) with U > is 
given by 

uv{x, q\y) = G(x, q\y)- ^^—7- r ^Oi, 9 1 2/), (3-22) 

l + yG(xi,g|a;i) 

where 

2 sin m sin ra/ 



G{x,q\y) = -> — - — . (3.23) 

7T z — ' + n z 

n=l 

Note that 

q 1 , v /"a/ 1 w 4 ^ sin(2n - l)y 

So [q y = / G(x, q\y)dx = - > — — . 

7 tt^ (2n - 1) + (2n - l) 2 ) 

According to equation (|3.20j) . the survival probability Sv(t \ y) is given by 



7T 







Sy{t I y) — I uv(x,t\y)dx (3.24) 
and the Laplace transform is 

§v(t I y) = I uv{x,t\y)dx. (3.25) 
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Using (J3.22)) . we find that the survival probabilities, without and with the 
Dirac killing, differ by 



So(q\y) ~ S v (q\y) 



VG(xx,q\y) 



1 + VG(xi, q | X\ 
To compute S (q \y), we use the formula 



-S {q I Xl)- 



(3.26) 



Su{q\z) = 



cos(nz) 



TV 



r* n z + q 

when x, y e]0,7r[. Then, 

G(x,q\y) 



cosh( v /gz) 1 
tanh( v /g7r) ^/qn J lyfq 

cos(y /= gz) 1 



tan( v /— g7r) yf—qit ) 2^— g' 



71 



for q > 
for g < 



Su(g|a; — y) — Su(q\x + y) 



7i 



Thus 



So(q\y) 



1 /cosh( v /g(x — y)) cosh( v /g(x + y)) 



„ 2^ V. tanh( v /gvr) 
1 



dx 



tanh( v /g7r) 
sinh( x /g(x - y)) - sinh( A /g(x + 



2gtanh( v /g7r) 

From a Taylor expansion around q = 0, we obtain that 

^o(g|y) = ^ + 0(^), 

where = — 37r 2 y — 7r3y 2 + y 3 , and similarly 

%|x 1 ) = ^ + 0( v ^). 



(3.27) 



(3.28) 



We conclude that in a bounded interval, the decay rate for the survival proba- 

67T 

bility of a free particle is exponential with a rate constant , which depends 

Q(y) 

on the initial position of the particle and is given by 

*(«I»)~«p{-^} fa «>!■ 

y(jr(Xl g ?/) 

For a given 1/ > 0, equation (J3.26j) contains the term - — , which 

1 + VG( Xl , q I Xl ) 
is 1 at q = 0, and at the first order approximation, when y 7^ x\, 

Sv(q I y) = S (q I y) - 5 (g I xi) = _ + O(^). (3.29) 
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We conclude that the strength V does not enter the first approximation of 
the survival probability, but the decay rate constant is bigger than in pure 
diffusion. More specifically, we obtain that 

S v (t | y) ~ exp ( - 6?i "^ 1 for t > 1. (3.30) 



The potential strength enters in the next term in the expansion of Sy. 

3.2.2 Computation of the conditional MFPT E[T \T <r,y] 

The conditional MFPT E[T \ T < r, y] is computed by using expression ()3.17j) 
as follows. Equation ()3.17|) . corresponding to the killing measure V5(x — X\), 
gives 



d 

E[T \ T < r,y] = - — ln{p(xi, q\y)} 



(3.31; 

g=0 



The Laplace transform of equation (j3.21|) with absorbing boundary conditions 
is given by 



2 V sm nx sm ny 

which gives for x — x\ 



u{x, q\y) = > ■ — ^ — g | y) + G(x, ? y, (3.32) 

7r ^y"' g + n^ 



u(xi, g|y) = — . , (3.33) 

I 21/ smnii smny 

7r f-' q + n 2 

and 

9 . . 9 , 9,1* 2^ sinnxi sinra/ 

— lnp(xi,g y) = F lnG ii,? j/ --r-ln H > ■ — 

oq oq oq \ 7i ^ q + n z 

= a(x 1 \y)+f3(x 1 \y) 



with 



a(zi|y) = ^lnG(xi,g| y) 



oo 

sm nxi sm 
ri A 

n=l 



oc 



5=0 sm nx\ sm ny 

2 



n=l 
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and 



. s d . ( „ 2V^ v~> sin nxi sin ny 



2K sin nxi sin ra/ 
7r ^ n 4 

n=l 

2 V sin n^i sin ra/ 

7T ' n 2 
1 

For x±, y G]0, tt[, it is well known that 

sin nx\ sin (7r — a?i)y 



2 



g=0 



7T z — ' n 2 7T 

71=1 

CO 



2 ^sinnxisinra/ x x y 2 2 2 . (^i+F 
— } - = [xl + y + 2tt ) 

7T ^— ' 7T 4 67T 6 



71=1 



7T 



so that finally, we obtain 

E[T\T<r,y] = -a(xi \y, ) + (3{x x \y) 

x\y{x\ + y 2 + 27r 2 ) — 7r(xf + y 3 ) 

6(7r — xi)y n + V(7r — xi)y 

3.3 Ratio measuring the distribution of particles 

According to the Fokker-Planck equation ([3.2)1 . the time dependent ratio i?(t) 
of the absorbed particles (particles leaving the domain, before being killed) to 
the killed particles at time t can be defined as 

/ J(x,t\y) ■ u{x) 

R(t) = ^ . (3.34) 

k(x)p(x, t\y)dx 



n 



More interestingly, we can define a steady state ratio Roo, which is the to- 
tal number of absorbed particles to the total number of killed particles after 
infinite time, by the expression 



J(x, 1 1 y) ■ v{x) dSx dt / J(x\y) ■ u(x) dSx 

' ' oc poo f r > \O.OOj 

\ \ k(x)p(x,t\y) dxdt / k(x)G(x\y) dx 

Jo Jn Jn 
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where G(x\y) is defined by the equation 

— p(y) — CG{x, | y) — k(x)G(x | y) for sjgU (3.36) 

with the forward operator C defined in eq. (|3.3|) . p(y) is the initial density, and 
J(x | y) is the flux density vector at point x, computed with respect to the 
function G(x | y). When p(y) = S(y), G is the standard Green function with 
boundary conditions given by equation ()3.7|) . 

We can define another ratio of interest: in a permanent regime, when 
a flux enters the domain through a part of the boundary, it is partitioned 
into the flux of absorbed and killed particles. When a steady state regime is 
achieved, we can define the ratio R s as above. We denote by dVti the part 
of the boundary, where a steady flux enters the domain. The steady state 
Fokker-Planck equation becomes 

= Cp(x | y) — k(x)p(x | y) for x,y G Q, (3.37) 

where the forward operator C is defined by (13. 3|) and the boundary conditions 



arc 



p(x \ y) = for , x G dVt, y G f2 a 
J(x | y) ■ v(x) = for x G Oil — dVt a - dQ h y G fi, t > 0. 

J{x | y) ■ v(x) = — $(x) for x G dVLi. 

The time independent flux is $(x) > 0. The external steady state flux of 
absorbed particles is 



J a = [ J{x\y)-u{x)dS x . (3.38) 



The total inward flux is 



J l= [ J(x\y)-u(x)dS x = [ §{x)dS x . (3.39) 
We define the ratio R s as 

/ J(x | y) ■ u(x) dS x / ®(x)dS x - k(x)p(x\y)dx 

R s = = z*k J J1 . (3.40) 

/ k(x)p(x \ y) dx i k(x)p(x \ y) dx 

Jn Jn 

The second part of the identity is a consequence of conservation of matter. 
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3.4 The one-dimensional case 



The fluxes Roo and R s can be explicitly evaluated in dimension 1, when the 
domain is a finite interval. The ratio R s was computed in ^H] i n the case of 
an interval [0,L], when the killing measure was uniformly distributed. 

We assume now that the inward flux at x = is a constant $ and at x = L 
an absorbing boundary condition is given. We consider here the case where 
the killing is a Dirac k(x) = k5(x — Xi), located at a single point x\ and k 
is a constant. The particles are only driven by diffusion, so the steady state 
equation ()3.37|) becomes 

D dp(x)_ _ k ^ p ^ = q for o < x < L 



dx 2 

dp(L) 



and the ratio is 



dx 

p(0) = 



d£(o) 



From an explicit computation of p(x), one can derive that 



fcp(xi) 



l + -{L-x x ) 



and 

^ = t./r D V ( 3 - 42 ) 

The result can be generalized to the case of a two hot spots in a straightforward 
manner. When killing occurs uniformly, the ratio R s decays as function of L 
like the function 1/ cosh(cL), (c = const.) This decay, compared to the 
decay of equation (|H.42jl . shows that any redistribution of the killing affects 
this ratio, which is discussed in the conclusion section. 

In the same spirit, we give an explicit expression for R^ in the case of a 
finite interval [0, L], where particles are free to leave the domain at the points 
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and L. Initially, we assume that the particles are located at a point x\. Here 
the killing occurs at the point y < x\. In that case the Green function G(x \ y), 
defined by ()3.36|) . is the solution of 

d 2 

-S(xi) = -D— - k(x)G(x\y) for x,ye[0,L] 

G(0\y) = 

G(L\y) = 
and by solving this equation, the ratio Roc is given by 
R = (L-x^k 



4 Conclusions, applications, and perspective 

We have provided in this paper a general mathematical framework to compute 
the distribution of "killed" and "absorbed" particles, after they flow into a 
bounded domain. The ratios Roo or R s of "killed" to "absorbed" particles 
are in general difficult quantities to estimate analytically However, in one 
dimension the exact dependency of the ratio on the geometry can be computed; 
we analyzed here two extreme distributions: a uniform distribution and a 
Dirac killing measure. Formulas (|3.42|) and 1/ cosh(cL), (c = const.) of [To] 
prove that the ratios depend on the killing distribution. For a general three- 
dimensional domain, R^ can only be estimated in asymptotic cases, where 
the absorbing boundary occupies a small portion of the boundary or when the 
support of the killing measure is small (see [T4"j). 

In the general context of microstructures in biological systems, the ratio R^ 
provides information about the total distribution of particles. When the killing 
measure is redistributed and a critical value of the ratio Roo is attained, new 
biophysical processes can be initiated that affect irreversibly the physiological 
properties of the microstructure. Indeed, if enough particles enter the structure 
and stay sufficiently long, they bind to a large number of molecules. When a 
critical number of bonds are made, a cascade of chemical reactions is initiated. 
Thus a threshold can be reached by simply redistributing the killing measure. 
The implementation of these changes at a molecular level is yet to be identified. 
The mean conditional time of being absorbed before killing, E(T,r < T), 
reveals not only the time spent inside the structure, but also how long it takes 
on the average for particles to arrive to a specific compartment. 

The spine-dendrite communication can be described in terms of quantities 
such as Roo and E(T,t < T). First, the regulation of calcium ions that 
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reach the dendrite can be achieved by various mechanisms. One possibility 
to decrease Roo is to increase the length of the neck, which really occurs in 
vitro experiments ^2]- I n that case, if the distribution of the killing measures 
is scaled with the dilation of the neck, the ratio changes, with no need 
to change the total killing measure (e.g., the number of pumps). A second 
possibility is to redistribute the killing measure in a way that affects the ratio 
Roo, as shown in our computations (e.g., from uniform to accumulation at a 
hot spot). We can predict from expression ([3.42)1 . that moving all the calcium 
pumps at the bottom of the spine neck reduces the number of ions arriving 
at the dendrite. Finally, the number of pumps can also be changed. All 
possibilities are expected to be observed and any particular choice should be 
understood in the context of its function. We expect that the distribution of 
pumps across the spine neck to be highly dynamic and driven by the mean 
electrical activity of the dendrite. In particular, we may wonder how such 
distribution changes in the wake of applying protocols such as LTP (Long Term 
Potentiation), which lead to long term changes at the synapse level JTU]- No 
results seem to be known about the effect of LTP on the pump redistribution 
in spines. In reality, as studied in the movement of ions inside the spine 
neck is not purely Brownian, but has a drift component, which affects the 
dynamics and changes the ratio. 

The mean time E(t | r > T) to arrive at the dendrite was used in JH] to 
confirm that calcium ions arriving at the dendrite originate at the spine head 
(not in external sources). This result is derived by comparing the experimental 
time scale with E(t \t > T). The mean time E(r | r > T) is thus a funda- 
mental parameter in the context of spine-dendrite communication, because it 
measures the mean time calcium ions enter the dendrite, and is related to the 
induction time of cascades of reactions, involved in modifying the synaptic 
weight. Changing the E(t | r > T) is a part of the spine regulation process. 
This can be achieved by various ways: changing the spine neck length, changing 
the number of pumps and their distribution. Various biological investigations 
(see for example jTT]) are dedicated to the elucidation of how such regulation 
is achieved at a biochemical level. 

Finally, the present computations assume that the neck width is small. If 
this is not the case, the one-dimensional approximation of the cylinder is no 
longer valid and pumps become insignificant. 
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